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FOREWORD 


The  program  described  here  consisted  of  determination  of  heating 
rates  and  diffusion  coefficients  for  the  turbulent  atmosphere  through 
photometric  and  photogrammetr i c  reduction  and  computer  analysis  of 
photographs  of  chemical  releases.  Related  earlier  work  by  the  Photo- 
Metrics  research  group  is  described  in  References  1,2,3  and  5. 

The  author  wishes  to  express  his  thanks  to  Mr.  Samuel  P.  Zimmerman 
(Contract  Manager)  of  AFGL  for  his  continued  encouragement  and  support, 
to  Z.  Reinhardt  and  J.  Lander  for  assistance  in  the  data  analysis,  and 
to  Dr.  l.L.  Kofsky  and  C.C.  Rice  who  also  made  important  contributions 
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SECTION  I 


INTRODUCTION 

The  objective  of  the  program  described  here  is  the  determination 
of  parameters  which  characterize  turbulence  in  the  upper  atmosphere. 
These  include,  but  are  not  limited  to,  turbulent  diffusion  coefficients 
and  atmospheric  heating  rates  derived  from  analyses  of  transport  and 
surface  brightness  distributions  of  tracer  chemical  clouds.  Micro- 
dens  i tometric  scanning  is  used  to  sample  the  radiance  distributions  of 
mainly  sunlight-scattering  artificial  clouds  which  have  been  recorded 
by  radiometrical ly  calibrated  photography.  These  brightness  data  are 
then  analyzed  by  digital  computer  methods  to  determine  both  Fourier- 
transform  and  configuration-space  statistical  descriptors  of  the  fluc¬ 
tuating  component  of  the  surface  brightness.  Features  of  the  one 
dimensional  energy  spectra  which  are  indicative  of  turbulence  are 
further  analyzed  to  determine  heating  rates  (rate  of  dissipation  of 
turbulent  kinetic  energy)  and  diffusion  coefficients  as  described  in 
Section  II.  Measurements  of  the  horizontal  wind  field  by  photographic 
triangulation  were  alsousedto  obtain  further  estimates  of  heating  and 
diffusion  (Section  III).  Support  tasks  and  an  inverse  triangulation 
scheme  are  discussed  in  Section  IV. 

A  summary  of  the  tracer  release  cloud  data  analyzed  is  shown  in 
T^ble  1,  and  a  partial  listing  of  results  from  the  "Aladdin  197V 
(through  the  night  release  series  --  Wallops  I.,  VA)  in  Table  2.  The 
results  of  these  and  other  analyses  were  reported  as  they  were  completed 
to  AFGL  scientists  for  their  immediate  use  in  the  interpretation  of 
atmospheric  turbulence  phenomena. 
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Table  1.  Summary  of  Analyses 
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SECTION  I  ' 


DETERMINATION  OF  HEATING  RATES  AND  EDDY  DIFFUSION 

Characterization  of  turbulent  flow  in  the  atmosphere  may  be  accom¬ 
plished  by  analyzing  surface  brightness  distributions  of  chen  i luminescent 
or  sunlight  scattering  tracers  released  from  rockets  (Pel  1-6).  Statis¬ 
tical  or  Fourier  methods  may  be  applied  to  the  mean  brightness  distri¬ 
bution  or  to  the  fluctuations  of  brightness  about  the  mean  distribution 
to  extract  information  about  the  structure  present. 

The  released  atoms,  molecules,  or  particulates  very  quickly  reach 
pressure  and  temperature  equilibrium  with  ambient  species  and  are  found 
from  exper:ence  to  conform  to  the  flow  field.  The  spatial  character¬ 
istics  of  the  flow  are  frozen  in  time  by  rad i omet r i ca 1 1 y  and  photo- 
grammetrical ly  calibrated  photography  from  ground  stations.  Analysis 
starts  with  microdensitometer  scans  across  the  images  of  this  frozen 
flow  field,  which  are  assumed  (by  the  ergodic  hypothesis)  statistically 
equivalent  to  the  time  resolved  output  of  a  fixed  physical  sensor  as  the 
flow  field  proceeds  across  it.  The  resulting  digital  density  informa¬ 
tion  is  converted  to  scene  brightness  through  the  HbD  characteristic  of 
the  film,  with  appropriate  corrections  arolied  ! o  account  for  reciprocity 
effects  in  exposure  duration,  spectral  sens  '  ■  >  •<  i  ty  of  the  film,  lens 
T-stop,  and  varying  sky  background  (Ref  3)  • 

Further  corrections  are  also  performed  to  account  for  the  transfer 
function  of  the  microdensitometer's  finite  sampling  aperture  and  for  the 
modulation  transfer  function  of  the  microdensitometer  itself,  which  is 
a  function  of  the  temporal  scanning  rate  of  the  instrument.  The  im¬ 
portance  of  utilizing  a  scanning  system  with  true  microdensitometer 
optics,  is  exp’ained  in  Ref  11.  Briefly,  to  obtain  the  precision  and 
accuracy  needed  for  these  analyses  the  response  of  the  instrument  must 
be  temporally  and  spatially  stationary.  Distortion  of  the  calculated 
power  spectra  can  result  from  nonuniform  reproduction  of  the  photo¬ 
graphically  recorded  dafa. 
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Photogrammet r i c  calibration  of  the  imagery  is  also  required  (next 
section)  to  determine  the  altitudes  of  the  points  of  measurement  and  to 
convert  film  plane  measurements  to  distance  scales  at  the  glow.  Time 
sequential  analyses  often  require  a  correction  of  scale  factor  due  to 
the  range  variation  produced  by  cloud  movement  by  the  mean  wind  field. 
Finally,  the  corrected  brightness  data  are  analyzed  by  Fourier  (or  other 
statistical)  methods.  Fourier  analysis,  or  calculation  of  one  dimensional 
turbulent  energy  spectra  is  carried  out  by  one  of  three  methods  -- 
standard  Fourier  analysis,  "fast"  Fourier  Transform  (FFT),  or  a  simple 
summation  process  (counting  method.  Ref  8)  which  determines  the  spectral 
energy  in  discrete,  nonuniform  wavenumber  bands.  Each  of  these  methods 
requires  large  amounts  of  input  data  to  achieve  a  desired  statistical 
accuracy  (Ref  10),  usually  20  to  50  parallel  microdensitometer  scans 
through  the  turbulent  region  of  interest.  Depending  upon  analysis 
method,  the  average  energy  spectrum  is  obtained  by  segmental  averaging 
or  by  averaging  of  individual  scan  spectra.  Examples  of  FFT  and  "count" 
energy  spectra  for  release  ETTY  are  shown  in  Figures  1  and  2. 

The  one  dimensional  scalar  energy  spectrum  E(k)  as  a  function  of 
wavenumber  k  may  be  utilized  to  determine  the  Fourier  characteristics  of 
turbulent  flow.  E(k)  often  exhibits  power  law  dependencies  on  wavenumber, 
the  most  well  known  being  the  Kolmogoroff  or  inertial  spectral  range 
where 

E  (k)  <r  k~5/3 

which  implies,  from  dimensional  arguments,  that  turbulent  energy  is  being 
fed  at  a  constant  rate  from  large  cells  into  small  ones  in  the  subrange 
of  cell  dimensions  in  which  this  slope  obtains.  A  further  application 
of  the  energy  spectrum  is  its  capability  of  isolating  features  of 
characteristic  wavelength  --  for  example  internal  atmospheric  gravity 
waves . 

When  an  inertial  spectral  range  is  present  the  heating  rate  or 
rate  of  viscous  dissipation  of  turbulent  kinetic  energy  e  may  be  found 


(1) 


A  L  3 

=  V3T7T7 

where  A  is  a  constant  E(k.)is  the  energy  at  wav -'number  !■  located  within 
:he  imrtial  subrange,  and  is  tne  scalar  :  jit  o'"  decay  o'"  .1-nsity 
fluctuations  (n")  due  to  nolecular  diffusion 
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0  =  the  molecular  diffusion  coefficient,  ard 

space  '•ate  of  change  of  turbulent  density 
fluctuations  of  tracer  material  in  the  release 


Jn 
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(— ■  )>  is  usually  determined  by  numerical  differentiation  of  the.  bright¬ 
ness  data  and  may  also  be  obtained  by  the  integral  formulation  (2). 
frequently  <  (4-  4  >  (and  thus  e)  Is  biased  toward  lower  values  because 

O  X 

of  limited  spectral  extent  in  the  integral  approximation,  or  by  inclusion 
of  smooth  regions  of  sky  background  in  the  numerical  method.  The  power 

law  form  of  E(k)  is  a  -5/3  to  0.2  k  where  a  transition  to  a  larger 

V 

negative  slope  occurs  (~3  or  greater);  k  -  \  is  the  limiting 

wavenumber  for  viscous  dissipation,  and  v  the  Kinematic  viscosity. 

Limbed  spatial  resolution  of  the  photography  determines  the  spectral 
extent  and  the  upper  integration  limit  of  equation  2.  If  this  limit 
is  less  than  0.2  ky  it  is  evident  that  <(3n'/dx  )>  may  be  substantially 
smaller  than  expected  and  estimates  of  e  will  also  be  low.  The  above 
method  has  been  used  successfully  to  determine  c  near  the  turbopuase 
(Ref  6  and  12).  If  the  energy  spectrum  extends  beyond  0.2  k  an 
accurate  determination  of  the  constant  A  in  (1)  may  be  made  (Ref  13). 

The  scalar  energy  spectrum  in  the.  inertial  range  is 


E  (k)  =  Ac 


-1/3  k  -5/3 


which  may  be  made  nond i mens iona 1  using  Kolmogorov  scaling 


1 1 


Ei(k)=-^  E(k)  =  k,  3  v  At'l/3  k'5/3  . 


v 


Further  scaling  the  spectrum  by  (■£-) 


5/3 


N(k)  =  Ak  4/3  vc‘1/3 

v 


or 


,3 , 
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“  u  A  3/1 

k  v  c 
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but 


-  k  “  v3 
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and 


A  =  N (k)  =  E(k) 


k  3v 
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( — — ) 
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Alternatively,  since  k  is  known,  and  assuming  the  constant  A  is  known, 
v  may  be  calculated  from 


v  =  . 


r  1/3 


(3) 


Equations  used  by  atmospheric  modelers  to  calculate  altitude  profiles 
of  v  become  increasingly  inaccurate  at  very  high  or  low  temperatures 
or  under  conditions  occurring  at  high  altitudes  (above  32  km),  standard 
atmosphere  calculations  of  kinematic  (and  dynamic)  viscosity  are  present 
tabulated  only  to  86  km.  Further  experimental  determination  of  the 
behavior  of  kinematic  viscosity  at  high  altitude  could  provide  the 
opportunity  to  refine  the  model  equations. 


Another  possible  formulation  analogous  to  the  spectral  form  for 
the  average  potential  temperature  may  be  used 


E  (k)  -aH  ^  c  2/3  k’5/3 
n  g  Dz 
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where  n  is  the  mean  density,  the  gradient  of  the  mean,  u  is  a  constant , 
and  g  the  gravitational  constant. 

One  further  method  of  calculating  t.  may  be  utilized:  The  average 

^2 

of  the.  square  of  the  turbulent  dt-’sity  fluctuations  n  is  calculated 
for  several  times  after  the  initial  release.  Then 

T 

_  _  3  (n"  ) 

n  St 


and  the  rate  of  turbulent  dissipation  is  calculated  from  equation  (1). 
This  formulation  for  decaying  turbulence  should  be  applicable  to  tracer 
pulse  or  trails  in  any  altitude  region. 

Evaluation  of  the  eddy  diffusion  coefficient  K  may  now  be  accom¬ 
plished  utilizing  turbulent  heating  rates.  Reference  7  shows  that  eddy 
diffusion  may  be  expressed  as 


K 


=  C  C 


1/3 


(5) 


where  C  is  a  constant  and  kQ  is  the  low  wavenumber  limiting  the  inertial 

spectral  range,  which  is  determined  by  inspection  of  the  turbulence 

spectrum.  Wavenumber  k  can  also  be  calculated  as 

o 

E(k)dk 

k  =  '  — - 

O  ^-,00 

I  k-1  E (k) dk 
'■'o 


if  the  spectral  range  experimentally  extends  beyond  0.2  k^. 

If  the  Brunt  Vaisala  frequency  N  of  the  atmosphere  is  known,  the 
eddy  diffusion  coefficient  may  be  calculated  as  (Ref  14) 

K  =  y  u/N2  (6) 

Alternatively,  Equation  (6)  can  also  be  used  to  determine  N*  when  K  and 
c  are  known. 
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Table  2.  Partial  Tabulation  of  Typical  Analysis  Results,  Alladin  197^ 


Richardson  number. 


ent  diffusion  coefficients  derived  from  the 
n  1 97^  series  of  releases. 


Figure  4.  Brunt  Vaisala  frequencies  calculated  for  the  Aladdin 
197^  series  of  releases. 


The  data  in  Table  2  (Alladin  197^'Wal lops  1.,  VA)  are  typical  of 
the  analysis  results;  l  was  determined  using  Equation  (1*)  and  K  from 
Equation  (5),  k  was  estimated  from  the  soectra  and  v  calculated  from 
Equation  (3)*  Figure  3  shows  the  calculated  eddy  diffusion  compared  to 

a  diffusion  model  with  turbopause  at  112  km  (Ref  15).  Figure  4  compares 

2  2 

determined  from  Equation  (6)  to  N  determined  from  1976  U.S.  Standard 

Atmosphere. 

Knowledge  of  the  atmospheric  wind  field  may  also  be  used  to 
determine  u  and  K  (Section  111),  or  time  sequences  of  low  frequency 
power  spectra  may  be  analyzed  by  methods  analogous  to  those  used  to 
determine  molecular  diffusion  (Ref  1),  however,  this  latter  method  is 
only  useful  when  the  tracer  cloud  retains  a  roughly  gaussian  structure. 
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SECTION  III 


PHOTOGRAPHIC  TR I ANGULAT ION 

The  use  of  photographic  photometry  to  determine  quantitative 
measures  of  atmospheric  turbulence  requires  that  accurate  photogrammetr i c 
calibration  be  performed  as  well.  Triangulation  determines  spatial 
locations  on  the  tracer  release  as  a  function  of  time  as  it  is  trans¬ 
ported  by  the  mean  wind  field.  The  position  of  the  glow  and  its  range 
determine  the  conversion  from  image  scale  to  dimensions  on  the  three- 
dimensional  object  itself.  Procedures  used  to  perform  these  calibrations 
are  detailed  in  References  16  and  17,  and  are  briefly  covered  below. 

Photogrammetric  triangulation  to  distant  objects  from  appropriately 
separated  ground  stations  is  based  upon  accurate  knowledge  of  each 
station  location  (in  latitude,  longitude  and  altitude  above  a  reference 
ellipsoid)  and  also  the  orientation  of  the  camera  film  plane  at  each 
station.  Determination  of  camera  orientation  is  accomplished  by  analysis 
of  photographs  of  star  fields.  Film  plane  coordinates  of  star  images 
and  equatorial  coordinates  of  the  stars  at  the  exact  time  of  exposure 
are  related  through  vector  equations  which  are  solved  to  determine  the 
exact  focal  length  of  the  camera  and  the  direction  of  the  camera  optic 
axis  in  right  ascension  and  declination,  and  also  the  rotation  of  the 
film  plane  needed  for  coincidence  of  the  stars  and  their  images.  The 
orientation  may  be  readily  expressed  in  any  other  coordinate  system  (e.g., 
horizon  or  geocentric)  by  rotations  about  the  coordinate  axes  which  effect 
the  coordinate  transformation. 

The  digital  film  plane  positions  of  the  tracer  cloud  as  viewed 
from  two  photographic  sites  are  then  related  to  determine  the  spatial 
location  of  the  glow.  Film  plane  coordinates  of  the  glow  are  determined 
using  an  automatic  scanning  densitometer  which  can  locate  and  digitize 
the  centerline  of  a  trail  or  puff  release  (Ref  18)  or  determine  the 
location  of  small  dots  drawn  on  the  negative  in  the  region  of  interest. 

The  latter  method  must  often  be  used  for  turbulent  portions  of  releases 
where  the  automatic  scanner  cannot  define  a  centerline  of  the  release. 
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Vector  lines  of  sight  for'  each  film  plan'  release  point  are  calculated 
and  an  iterative  procedure  determines  the  best  match  of  each  Site  1 
point  with  a  Site  2  point.  The  spatial  position  of  the  trail  is  deter¬ 
mined  by  calculating  the  latitude,  longitude  and  altitude  of  the  inter¬ 
section  of  the  vector  iint-s  of  sight  of  matched  points.  Some  of  the 
problems  which  arc  associated  with  the  i mp I einen t a t i on  of  the  method  are 
detailed  in  Reference  19.  An  inverse  of  thi'  t r i angu I  a t i on  method  was 
used  to  approximately  project  altitude  grids  onto  the  film  plane  records 
of  older  releases  given  the  known  wind  profile  (Section  IVj. 

Measu remen t s  of  the  horizontal  wind  field  are  accomplished  by 
analysis  of  the  motion  of  the  trail  for  several  frames  over  total 
durations  of  one  to  two  minutes.  Least  squares  analysis  of  trail  position 
as  a  function  of  time  at  each  altitude  yields  the  average  positional 
change  and  average  velocity  components. 

Further  detail  of  atmospheric  turbulence,  heating  rates,  and 
diffusion  may  be  studied  with  the  aid  of  these  horizontal  winds. 
Richardson  number  shown  in  Table  2,  which  is  a  measure  of  the  local 
stability  or  instability  of  the  atmosphere,  was  determined  from 

R.  =  n2/ (;*v/nz) 2 

t 

where  3V/3Z  is  the  vertical  shear  of  the  horizontal  wind.  Usually 
R  }  indicates  stability,  R.  •  instability  or  turbulence.  Whenever 
winds  are  used  as  a  basis  for  measuring  turbulent  heating  and  diffusion 
the  value  of  R.  must  be  less  than  i  in  the  region  of  the  measurements. 

The  shear  used  for  the  data  in  Table  2  is  an  average  over  two  km  alti¬ 
tude  segments;  a  smaller  increment  in  alt'tude  would  produce  a  larger 

shear  and  decreased  R..  The  energy  deposition  rate  per  unit  mass  due  to 

1  2 

the  vertical  shear  of  the  horizontal  wind ,  c-  ?  v( 3v/3z)  is  an  additional 
measure  of  heating.  For  turbulent  layers,  turbulent  velocity  fluctuation 
v'  about  the  mean  velocity  may  be  calculated  and  heating  may  be  expressed 


l  =  A'"  V°  -/f 


where 

order 

equat 


where 


is  the  outer  length  of  the  turbulence  and  A'  is  a  constant  of 
unity.  Under  certain  conditions  the  one  dimensional  energy  balance 
on  may  be  solved  to  provide  a  vertical  diffusivity 


K(z) 


c 

Ov/3z)2 


a  K 

m 


K  is  the  momentum  transfer  coefficient, 
m 
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SECTION  IV 


SUPPORT  TASKS 

The  work  described  in  this  Section  relates  direct  1/  to  determining 
reliable  turbulence  spectra  at  a  series  of  known  atmospheric  altitudes. 

It  involves  choosing  photographic  records  of  past  releases  which  l)  show 
turbulent  structure,  and  2)  have  the  appropriate  photogrammet r i c  and 
radiometric  calibrations,  or  are  amenable  to  reconstruction  of  these 
calibrations,  for  quant i tat i ve  analysis. 

First,  we  surveyed  and  cataloged  a  library  of  film  records  of 
AFGL  chemical  releases  dating  from  1962  to  197*i.  Of  an  estimated  1/00 
rolls  of  film  that  were  available  for  this  purpose  866  received  pre¬ 
liminary  examination  to  determine  those  releases  which  could  be  candi¬ 
dates  for  useful  analyses.  Of  these  235  were  judged  to  contain  no  evi¬ 
dence  of  atmospheric  turbulence,  and  another  319  were  found  to  be  dupl i 
cates  (more  than  one  camera  is  usually  operated  at  each  photographic  site). 
The  remaining  312  films  were  carefully  examined  and  numbered,  and  a  loq 
of  data  prepared  containing  the  following  information:  Log  Number, 

Release  Name  or  identification,  Date,  Site.  Camera  Number,  Film  Type, 
Wavelength  Filter,  Launch  Time,  Release  Material,  Type  of  Release  (trail 
or  puff),  whether  a  calibration  stepwedge  had  been  made,  whether  the 
film  had  been  previously  used  for  t r i angulat ion ,  frame  numbers  or  time 
code  where  turbulence  was  present,  and  comments  or  description  of  the 
turbulence. 

Of  the  films  in  the  catalog  77  were  found  which  had  turbulence 
present  and  also  had  the  necessary  calibrations  for  meaningful  analyses. 

The  amount  of  data  which  was  suitable  for  analysis  was  considerably  more 
then  could  be  analyzed  during  the  course  of  the  work  and  the  cataloging 
procedure  was  halted  at  this  point.  Photographic  prints  of  the  turbulent 
areas  were  made  which  could  be  annotated  to  preserve  a  record  of  areas 
which  were  digitized  for  analysis;  as  this  is  part  of  the  working  data, 
these  images  are  not  reproduced  here. 
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A  major  problem  was  encountered  with  these  older  data  (all  pre-1975). 
All  of  these  releases  chosen  for  analysis  had  been  previously  triangulated 
to  determine  release  position  as  a  function  of  time  and  altitude,  and 
wind  profiles  had  been  calculated  from  these  positions.  It  was  expected 
that  we  would  correlate  the  computer  printouts  of  trail  position  with 
the  t r i angu 1  at i on  photographs,  on  which  were  marked  a  1 1 i tude- ident i f y ing 
dots,  to  determine  the  altitudes  of  regions  to  be  analyzed.  However, 
these  printouts  proved  to  be  no  longer  available  either  from  AFGL  or 
from  the  contractor  who  had  performed  the  wind  determination.  The  only 
printout  data  available  were  tabulations  of  horizontal  wind  profiles 
contained  in  reports.  A  lengthy  procedure  (a  sort  of  inverse  triangula¬ 
tion)  was  programmed  to  circumvent  this  problem  by  utilizing  known  or 
estimated  camera  site  and  orientation  data,  and  estimated  release  rocket 
trajectories. 

Tabular  wind  profiles  were  first  entered  into  computer  files  at 
one  kilometer  altitude  intervals.  Then  using  the  known  site  location, 
approximate  camera  viewing  angles,  and  an  approximate  latitude,  longitude 
location  for  the  release,  the  computer  program  created  a  model  release 
along  an  interpolated  rocket  trajectory  (from  the  few  position  data  points 
available),  and  then  transported  this  release  volume  using  the  known 
wind  profile.  At  predetermined  increments  of  time  after  injection  of  the 
tracer,  the  program  utilized  vector  t r i angu 1  at  ion  methods  to  reconstruct 
and  plot  at  1  km  intervals  the  image  plane  view  for  the  particular  camera 
used.  The  plotted  view  could  then  be  overlaid  on  the  film  to  determine 
a  1 1  i  tudes . 

The  estimated  error  for  this  procedure  proved  to  be  large  since 
some  of  the  input  parameters  must  be  estimated.  If  the  site  location, 
camera  orientation,  and  rocket  trajectory  were  all  known,  the  error 
would  be  on  the  order  of  )  km  in  altitude.  The  resulting  error  for  the 
data  analyzed,  however,  appears  to  be  closer  to  2  to  1*  km. 
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SECTION  V 


CONCLUSIONS 

Calculation  of  turbulent  atmospheric  heating  rate:  and  turbulent 
diffusion  coefficients  were  made  from  data  extracted  from  photographic 
records  of  releases  in  the  stratosphere,  mesosphere,  and  lower  thermo¬ 
sphere.  Analysis  methods  which  proved  most  useful  were  Fourier  energy 
spectra  to  characterize  this  turbulence,  numerical  analysis  to  determine 
the  gradient  of  the  mean,  and  analysis  of  horizontal  winds  and  shears. 
The  combination  of  spectrally  derived  parameters  coupled  with  the 
measurement  of  the  gradient  of  the  mean  has  produced  the  most  reliable 
determinat ions  of  heating  and  diffusivity.  Improvements  in  spectral 
technique  have  allowed  the  spectral  regime  transition  at  0.2  k^  to  be 
determined  and  thus  v  may  also  be  found.  The  Brunt-Va i sa 1  a  frequency 
may  also  be  evaluated  directly  from  the  spectral  data  without  recourse 
to  temperature  measurements. 

Calibration  data  for  many  of  the  older  releases  (pre  1975) 
were  found  to  be  only  marginally  adequate,  as  reconstruction  of  the 
calibrations  was  time  consuming  and  also  subject  to  errors  on  the  order 
of  2  to  b  km  in  altitude.  A’though  this  was  a  disappointing  finding, 
approximately  50  percent  of  these  trails  do  ha\  ;  sufficient  calibration 
information  to  determine  altitudes  to  approx  lately  one  km  accuracy. 

Results  of  data  analysis  reported  here  (refer  to  Tables  1  and  2) 
are  representative  only,  and  do  not  reflect  the  large  volume  of  data 
reduced  and  analyzed,  whose  results  were  reported  to  AFGL  scientists  as 
they  were  obtained. 
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